Structural properties of dense hard sphere packings 
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The structural properties of dense random packings of identical hard spheres (HS) are investigated. 
The bond order parameter method is used to obtain detailed information on the local structural 
properties of the system for different packing fractions <j>, in the range between <j) — 0.53 and 
cj) — 0.72. A new order parameter, based on the cumulative properties of spheres distribution 
over the rotational invariant wq, is proposed to characterize crystallization of randomly packed HS 
systems. It is shown that an increase in the packing fraction of the crystallized HS system first 
results in the transformation of the individual crystalline clusters into the global three-dimensional 
crystalline structure, which, upon further densiflcation, transforms into alternating planar layers 
formed by different lattice types. 



The model of hard spheres (HS) is of fundamental im- 
portance in condensed matter physics and material sci- 
ence since it successfully reproduces the essential struc- 
tural properties of liquids, crystals, glasses, colloidal sus- 
pensions and granular media. Packings of HS have also 
been used in solving important problems of information 
and optimization theories [HQ. In this Letter we fo- 
cus on the structural properties of dense three dimen- 
sional (3D) HS systems at different packing fractions 
4> = f per 3 , where p = N/V is the density of N hard 
spheres in a system volume V and a is the diameter of 
the spheres. In particular, we take a large set of pack- 
ings composed of TV = 10 4 identical spheres with periodic 
boundary conditions, generated using Jodrey-Tory (JT) 
[H|, 0] and Lubashevsky-Stillinger (LS) [3] algorithms as 
described in detail in Refs. 0, H[. The corresponding 
packing fractions vary in the ranges <p ~ 0.53 — 0.71 
and (j) = 0.58 — 0.68 for JT and LS protocols, respec- 
tively. Both ranges include the random close packing 
state (RCP) at </> c ~ 0.64 (Bernal limit) @. The main 
purpose is to look into the details on how densiflcation 
of the disordered solid state affects the local and global 
order of the HS system. 

To define the local structural properties of the sys- 
tem we use the bond order parameter method [9], which 
has been widely used in the context of condensed matter 
physics [9, 10], HS systems HM20L complex plasmas 12 ll - 
[25j, colloidal suspensions [26l430j , granular media [3 if . 
etc. In this method the rotational invariants of rank I of 
both second qi(i) and third wi(i) order are calculated for 
each sphere i in the system from the vectors (bonds) con- 
necting its center with the centers of its N nn (i) nearest 
neighboring spheres: 
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Yj are vectors connect- 
ing centers of spheres i and 7. In Eq.© 

" [mi m 2 ms 

are the Wigner 3j-symbols, and the summation in the 
latter expression is performed over all the indexes rrii = 
satisfying the condition mi + m 2 + ms = 0. In 
3D the densest possible packings of identical hard spheres 
are known to be face centered cubic (fee) and hexagonal 
close-packed (hep) lattices (<j) ~ 0.74). To detect these 
structures we calculate the rotational invariants #4, q$, 
and wq for each sphere using the fixed number of nearest 
neighbors N nn = 12. A sphere whose coordinates on the 
plane (#4, (fa) are sufficiently close to those of the ideal 
fee (hep) lattice is counted as fcc-like (hep-like) sphere. 

It has been shown recently [32[ that the cumulative 
distribution function of the form 
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is an extremely sensitive measure of the structural or- 
der in the system. Here u(wq) is the distribution of 
spheres over the rotational invariant wq normalized to 
unity. The relevant order parameter is the position of the 
half- height of the cumulative distribution WQ h , so that 
We(wQ h ) = 1/2. Let us check whether this order param- 
eter can be used to identify the onset of crystallization 
in the randomly packed HS system. 

Figure [1] shows the dependence of the order param- 
eter WQ h and the mean value of q^ [defined as (q^) = 

N 

(l/N) #6 (&)] on the packing fraction <j>. Both proto- 

i=l 

cols used here (JT and LS) reveal nearly identical struc- 
tural properties as clearly seen from Fig. [TJ Explosive- 
(1) like growth of the order parameter WQ h at (j) > (j) c re- 
flects the appearance of crystalline order (hep-like and 
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FIG. 1: (Color) Order parameter Wq h (red and blue lines cor- 
respond to JT and LS packing protocols, respectively) and 
the mean value of the rotational invariant (qe) (green) ver- 
sus packing fraction cj) for the dense randomly packed HS 
system. Explosive- like growth of the parameter w^ h ((j)) at 
cj) > cj) c ~ 0.64 is a signature of the appearance of crystalline 
order (hep-like and fcc-like spheres) in the system. On the 
other hand, increase in (qo) near cj) c is monotonous and slow. 
Inset (a) shows distributions of spheres over the parameter 
wq for the two packing fractions (<fi ~ 0.655, blue solid curve; 
and <j> ~ 0.645, red solid curve). Cumulants of these dis- 
tributions are plotted by dashed curves of the corresponding 
color. Emergence of the peak in the distribution over wq cor- 
responds to the onset of crystallization in the system. The 
values of wq for perfect hep and fee lattices are also indicated 
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Insets (b, 



c) demonstrate two realizations of sphere packing with very 
close packing fractions near cj) « 0.7, which have quite differ- 
ent {qo) values. Example of disordered packing at cj) ~ 0.55 is 
shown in the inset (d). The spheres are color-coded by their 
qe value (red and green colors correspond to fee- and hep- 
like spheres, respectively), revealing fec-dominated (b) and 
hep-dominated (c) packings. The possibility for the dense HS 
systems to have quite different composition (with different 
relative density of fee and hep crystalline configurations) is 
the main reason behind the pronounceable oscillations in (qe) 
observed for 6 > 0.68. 



fcc-like spheres). This proves the extreme usefulness of 
this structural measure to locate the onset of crystalliza- 
tion, as demonstrated here in the context of randomly 
packed HS systems. Inset (a) shows probability distri- 
bution functions (PDFs) of wq for the two packing frac- 
tions [(j) c± 0.645 (red curve) and <\> c± 0.655 (blue curve)]. 
The peak on the presented PDF corresponding to the 
denser system (blue curve) corresponds to the appear- 
ance of hep-like spheres. Cumulative distributions for 
these PDFs are also plotted making clear the basic idea: 
It is very convenient to employ smooth cumulative dis- 
tributions of this kind to characterize structural prop- 
erties of systems near crystallization. The dependence 
of the mean q§ value on <fi shown in Fig. [1] also demon- 




FIG. 2: (Color). Probability distribution functions (PDFs) 
of the rotational invariant q§ for different packing fractions 
<j> in the regimes corresponding to the onset of crystalliza- 
tion (a) and to the dense crystallized HS systems (b). Top 
panel (a) shows the PDFs for ~ 0.645 (blue), ~ 0.65 
(cyan), <f> c± 0.655 (green) and <j> c± 0.66 (red). The val- 
ues of qe for perfect hep and fee lattices are also indicated 
(<?6 Cp ~ 0.48, q { 6 cc ~ 0.57). Bottom panel (b) presents the 
PDFs for cj) ~ 0.685 (red), ~ 0.695 (green), ~ 0.71 
(cyan), ~ 0.715 (blue) and cj) c± 0.716 (magenta). In the 
latter case quasi-crystalline (QC) spheres with q^ c ~ 0.52 
(qQ Cp < q^° < ql cc ) can be identified (the typical QC ar- 
rangement is shown in the bottom panel; the cluster looks 
like torsional modification of hep/fee lattices 22]). The cu- 
mulants Qq of these distributions, Qq(x) = n(qe)dqe are 
plotted by dashed lines of the same colors. These curves, 
in particular, demonstrate abundance of the hep- and fcc-like 
spheres and reveal the oscillatory behavior of nf cc and nh cp for 
sufficiently dense packing (cf. variations of (#6) for cj) > 0.68 
in Fig. [J). 



strates interesting properties. Relatively slow growth of 
(qe) for cj) < <j) c is associated with the appearance of 
single hep-like spheres [19[. For </> > </> c the parameter 
(qo) grows more rapidly. Further densification of the sys- 
tem results in strong oscillations of which start at 
<p ~ 0.68. Geometry gives a simple explanation of such 
an oscillatory behavior: The point is that the hep and 
fee crystal structures are characterized by quite different 
values of the parameter q$ (for ideal hep and fee lattices 
gg Cp « 0.48 and q^ c « 0.57). Dense HS packings can 
be apparently realized (for a given </>) with significantly 
different relative quantity of hep and fee phases [see in- 
sets (b) and (c) in Fig. [TJ demonstrating hep-dominated 
and fec-dominated packings, respectively). This explains 
the oscillations. Note that similar results were obtained 
recently in 1<J [20[ . Returning to the regime of packing 
fractions near C , Fig. [1] clearly demonstrate that the or- 
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FIG. 3: (Color). Relative density of hep-like and fcc-like 
spheres, nh cp (green) and nf cc (red), as well as the total 
density of crystalline spheres n to t = ^hc P + nf cc (blue) , ver- 
sus the packing fraction (j). For moderate packing fractions 
((/) < 0.645), the increase of ntot with <j) is due to the appear- 
ance of single hep-like spheres. The increase has a power- 
law character (ntot oc a , where a c± 10). Much more 
rapid power-law growth of both n to t and nf cc takes place for 
(f) > 0.645 (ntot oc <j>P , nf cc oc 7 , where /3 c± 50 and 7 ~ 100). 
Transition from hep- to fec-dominated packing (nh cp < nf cc ) 
occurs at c± 0.66. Oscillations of nf cc and nh cp which are 
clearly seen for sufficiently dense systems (<fi > 0.68) reflect 
two preferable possibilities of HS packing (cf. Figs. [T] and [2]). 
The dependence of the order parameter WQ h on <j> is also plot- 
ted (dashed black line) to indicate the onset of crystallization. 



der parameter Wq is considerably more sensitive to the 
onset of crystallization in the HS system as compared to 
the mean value of the rotational invariant (qe). 

Figure [2] shows the probability distribution function of 
the order parameter q$ for different packing fractions cor- 
responding to the onset of crystallization (top panel) and 
dense crystallized packings (bottom panel). Cumulative 
distributions Qq(x) = ^ Tt{qo)dqQ are also shown for 
each PDF. The relative numbers of spheres in hep- like 
(nhcp) and fcc-like (nf cc ) state are plotted in Fig. [3] as 
functions of <p. The total number of spheres involved 
in crystalline arrangements (n to t = ^hep + ^fec) is also 
shown. For <\> < <p c the relatively slow increase of n to t 
with (j) is associated with emergence of single hep-like 
spheres, which are nearly randomly distribute over the 
system volume. The increase has a power-like character 
(n to t oc (j) a with a ~ 10). At <p ~ <p c the slopes of n to t(0) 
and nf cc (0) drastically increase. The number of fcc-like 
spheres reaches that of the hep-like spheres (nf cc ~ nhcp) 
at <j) ~ 0.66, where transition from hep-dominated to fec- 
dominated packing occurs. Oscillations of nf cc and nhcp, 
clearly seen for <j> > 0.68 correspond to the property of 
dense HS packings discussed above (see Figs. Q] and [2]) 
and also observed in Ref. [19] with the help of Delaunay 
simplexes. 
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FIG. 4: (Color online). Relative fraction of fcc-like spheres, 
^fcc/(^fcc + ^hep), for dense HS packings (0 > 0.68). About 
35 realizations generated by the JT algorithm are analyzed. 
Two peaks reflect domination of either fee or hep phases in 
the system. 



Figure |4] shows the probability distribution of the ra- 
tio nf cc /(nf cc + nhcp) for the dense (</> > 0.68) HS sys- 
tems. Two peaks of the distribution indicate that dense 
HS systems prefer to be dominated by either fee phase 
(with somewhat higher probability) or by hep phase 
(with somewhat lower probability). Equal distributions 
between fee and hep phases are considerably less proba- 
ble. However, the lack of statistics (only 35 packings were 
available for the analysis) and specifics of the JT algo- 
rithm may also affect this observation. For this reason 
we do not elaborate further on this point here. 

Cluster analysis of the arrangements of hep-like and 
fcc-like spheres reveals another interesting structural 
property of dense HS packings. Figure [5] shows how the 
distribution (and the shape) of clusters composed of fcc- 
like spheres varies upon densification of the system. For 
relatively low volume fractions, <j) < 0.65, fcc-like clusters 
consist of only few spheres. Increase of (j) results in the 
appearance of several big (containing few tens of spheres) 
fcc-like clusters. A typical cluster has a complicated 3D 
shape (see Fig. [5] for <p ~ 0.66). Eventually, nuclei of fcc- 
like and hep-like clusters transform into global crystalline 
structure occupying most of the system volume. Further 
increase in (f) leads to a structural change: a transition 
from essentially 3D character of this structure to 2D lay- 
ers composed of fcc-like and hep-like spheres occurs. This 
transition is illustrated in Fig. [5] (for <fi ~ 0.68). An ex- 
ample of such "sandwich-like" structure can be also seen 
in Fig. OJc). 

To conclude, we have demonstrated that the order pa- 
rameter, based on the cumulative distribution of spheres 
over the values of the rotational invariant wq, is a very 
sensitive measure of the onset of crystallization in the 
system of randomly packed hard spheres. The proposed 
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FIG. 5: (Color). Distribution of fcc-like arrangements for 
four different packing fractions <j> of the dense HS system. 
Spheres are color-coded by the value of the rotational invari- 
ant qe (red color corresponds to fcc-like spheres). At low 
(j> only a few individual fcc-like spheres and small clusters 
from fcc-like spheres can be seen. Densiflcation leads to the 
formation of large 3D clusters of fcc-like spheres. Eventu- 
ally, clusters of fcc-like and hep- like spheres transform into a 
global crystalline structure. Further densiflcation results in 
the structural transition, when essentially 3D shape of this 
global crystalline structure suddenly changes to planar 2D 
layers formed by fcc-like (shown in the figure) and hep-like 
spheres. 

order parameter is especially convenient indicator of the 
appearance of hep and fee crystallites in the system. 
We used this to investigate how the relative distribution 
between hep and fee arrangements varies with increas- 
ing the sphere packing fraction. We observed that at 
(j> ~ 0.68 structural transition occurs: The essentially 
three dimensional shape of the global crystalline aggre- 
gate suddenly changes to planar sandwich-like (layered) 
structure. 
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